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ABSTRACT 

We present new cloudy and cloudless model atmospheres for brown dwarfs 
using recent ab initio calculations of the line list of ammonia (NH 3 ) and of the 
collision-induced absorption of molecular hydrogen (H2). We compare the new 
synthetic spectra with models based on an earlier description of the H2 and 
NH 3 opacities. We find a significant improvement in fitting the nearly complete 
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spectral energy distribution of the T7p dwarf Gliese 570D and in near infrared 
color-magnitude diagrams of field brown dwarfs. We apply these new models to 
the identification of NH 3 absorption in the H band peak of very late T dwarfs and 
the new Y dwarfs and discuss the observed trend in the NH 3 — H spectral index. 
The new NH3 line list also allows a detailed study of the medium resolution 
spectrum of the T9/T10 dwarf UGPS J072227.51-054031.2 where we identify 
several specific features caused by NH 3 . 

Subject headings: stars: low mass, brown dwarfs — stars: atmospheres 



Introduction 



Brown dwarfs, with masses between ~ 0.01 - 0.075 M , are not massive enough to 
maintain steady hydrogen fusion in their cores and, except when they are very young, have 
effective temperatures below that of main sequence stars. Their low atmospheric tempera- 
tures allow the formation of molecules and condensates that leave conspicuous imprints on 
their spectral energy distribution (SEP ) that have led to the creation of two spectral classes: 
the L dwarfs with T e g ~ 2400 - 1400 K (jKirkpatrickl 120051) and the cooler T dwarfs extending 
down to T e g ~ 500 K (ILeggett et al.ll2010al ). The recent discovery of several extremely cool 
dwarfs by the Wide-field Infrared Survey Expl orer (WISE) motivat es the creation of a new 



spectral class: the Y dwarfs with T e g < 500 K (jCushing et al.ll201ll ). 



The spectra of L dwarfs are characterized by absorption bands of V O, TiO, H?Q, CO, 



FeH, CrH and the very broad doublets of Na I (0.59 fim) and K I (0.77 fim) (IKirkpatrick et al. 



19991 ). While condensates, primarily corundum, magnesium silicates and iron, do not have 
spectral features as distinctive as those of molecules, their continuum opacity reddens the 
near-infrared spectrum of L dwarfs significantly. The effects of the cloud opacity peak 
at T c fj ~ 1500 -1600 K or spectral type L6-L8. In the cooler T dwarfs, the condensate 
cloud layer si nks below the photosphere, TiO and VO have been lost to Ti0 2 , V0 2 and 
condensates ( Burrows fc Sharp! 19991: Loddersl 2002 ) the carbon chemistry has shifted to 



favor CH 4 (IBurgasser et al 



mid-infrared around T2 (jCushing et al 



2002; 



Lodders fc Feglevl 120021 ) , and NH 3 bands appear in the 



20061 ). In the ve ry latest T dwarfs and in the new Y 



dwar fs, weak bands of NH 3 appear in the near-infrared (IDelorme et al.ll2008l ; ICushing et al. 
201lh . 



Due to their additional degrees of freedom and greater complexity, the opacities of 
molecules are not nearly as well understoo d as those of atoms and ions. Th e primary molec- 
ular opacity data bases are the HITRAN jRothman et alil2003l . I2OO5L I2OO9I ) and the GEISA 
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( IJacquinet-Husson et al.l Il999l . 120051 ) compilations, which were developed for modeling the 
Earth's atmosphere and consider the opacity of gases at temperatures of 300 K and below. 
Since the temperature in brown dwarf atmospheres can be up to 10 times higher, those line 
lists are missing the so-called "hot bands" that arise from the absorption of a photon from 
an excited lower level that is not thermally populated at 300 K. In some cases, the molecular 
line spectrum is so rich that individual transitions are unresolved in laboratory spectra or 
it becomes impossible to assign the observed transitions to specific molecular energy lev- 
els. The line lists are thus incomplete for the conditions encountered in brown dwarf atmo- 
spheres, especially for higher photon energies, hot bands, and weak transitions corresponding 
to near infrared and optical wavelengths. Incomplete molecular line lists li mit the predic- 



tive capability of atmosphe re models and synthetic spectra of brown dwarfs (jCushing et al. 



20081 ; IStephens et al.l 120091 ) , a problem that was recognized in the first attempts to model 



the observed spec trum of a brown dwarf, the T7p dwarf Gliese 229B (IMarley et al.l 11996 



Allard et al.lll9961 ). In view of the difficulties in analyzing laboratory spectra of high temper- 
ature molecular gases, an alternative is to compute a molecular line list with a first-principles 
approach. 

The past 15 years of brown dwarf research have motivated the theoretical modeling of 
the opacity of several key mole cules. New molec ular li ne lists for HqO dPartridge Schwenke 



1997 



Barber et alJboOffl . CO JGoorvitchlll994h. TiO JSchwenke!ll998[ ). VP Jpiezlll998h . FeH 



Dulick et al.l 120031 ) , CrH ( iBurrows et al.l 120021 ) , have improved the quality of the brown 
dwarf model spectra significantly. Methane (CH 4 ) has been and remains the most important 
molecule obs erved in brown dwarf spectra (defining the T spectral class, but appearing as 
early as L5; ( jNoll et al.ll2000l )) with a ver y incomplete line list at near infrared w avelengths 
although some progress has been made (IFreedman et al.l 120081 ; ISchwenkd |2002| ) . The two 
most recent developments in this area are new calculations of a line list of ammonia (NH 3 ) 
and of the continuum opacity from collision-induced absorption (CIA) by H 2 molecules. 

We study the effect of these two much anticipated new opacity calculations on the 
spectra and colors of brown dwarfs. In section §2, we discuss the nature of H2 CIA opacity, 
how the new calculation changes the opacity and how if affects the SED of cool brown dwarfs. 
The role of the NH 3 opacity in brown dwarf spectra is presented in section §3, where we 
show how the new line list for NH 3 affects primarily the near-infrared spectrum and brings 
new spectral features that can help define the T/Y dwarf boundary. We demonstrate the 
resulting improvement in synthetic spectra by comparing with the well-characterized SED 
of the T7.5 dwarf Gliese 570D and near infrared color- magnitude diagrams in section §4. 
The improved NH 3 line list provides the opportunity to comment on the role of ammonia in 
extremely cool dwarfs near the T/Y transition and on recent detections of this molecule in 
the spectra of extreme objects (section §5). Concluding remarks are provided in section §6. 
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2. Collision-induced absorption by H2 

Molecular hydrogen H2 is the most abundant component of brown dwarf atmospheres. 
Because of its high degree of symmetry, the H2 molecule has no permanent dipole moment 
and only quadrupole transitions are allowed, which are very weak. During a collision with 
another particle, however, the compound H 2 -particle system has a temporary dipole moment 
and can thus undergo a radiative dipole transition. The collision-induced absorption pro- 
cess is significant only at relatively high densities where collisions are frequent, such as in 
the atmospheres of planets and brown dwarfs. Three-body and higher order collisions also 
contr ibute to CIA but are probably negligible in brown dwarf atmospheres (IHare fc Welsh 
1958h . 



While the probability of a dipole transition in a typical molecule involves only molecular 
properties, CIA also depends on the frequency of collisions with other atoms or molecules 
(primarily other H 2 and He, by far the most abundant species), which is linearly dependent 
on the density of the gas. For a given metallicity, the relative importance of CIA absorption 
in a brown dwarf atmosphere is a measure of the density of the gas, which is closely tied 
to the surface gravity. On the other hand, the importance of H2 CIA decreases in higher 
metallicity atmospheres as the absorption from molecules involving heavier elements (CH4, 
CO, H 2 0, NH3) becomes dominant. This sensitivity of Ho CIA to gravity and metallici ty is 
a useful spectral diagnostic in T dwarfs (IBurgasser et al.l 12002 . l2003al ; iBurgasserl l2004bl ) . 



The H 2 CIA transitions occur between vibrational-rotational levels of the H 2 molecule 
and form a series of bands corresponding to At; = 0, 1, 2, 3, etc., where v is the vi- 
brational quantum number of the molecule. Each band is formed by transitions between 
the different rotational levels for a given Av. Collisional broadening of the individual 
lines is comparable to the line spacing and thus H 2 CIA is effectively a continuum source 
of opacity. Borysow and her collaborators calculated the H 2 -H 2 and H 2 -He CIA coeffi- 
cients in the temperature regime of interest for brown dwarf and white dw a rf atmospheres 



(Bor vsow fc Frommhold 



2001 



Borysow 



200 



2; 



1986 



Borvsow et al. 



1989 



Gustafsson fc Frommhold 



brown dwarf modelers (IFreedman et al 



200 



.Jorgensen et al.l 120001 ; iBorysow et al. 
and those are universally used b y 



Sharp fc Burrows! 120071 : lAllard et al.ll2001[ ). 



The H 2 -H 2 and H 2 -He CIA coefficients are shown in Figures 1 and 2, respectively. While 
these two sources of opacity have similar frequency (u) and temperature (T) dependences, 
H 2 -He plays a secondary role since the collisions with He are about 10 times less frequent 
than with H 2 as follows from their relative number densities. The absorption coefficient per 



1 Additional H? CIA references as 



well as codes and tables can be found at 



http : / / www . astro . ku . dk/ ^aborysow/programs/index . html 
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H 2 molecule is 

k H2 (u,T) = n H2 a H2 _ H2 (v,T) + n nc cr H2 _ He (v,T), (1) 

where is the number density of species i and cr H2 _i is the H 2 — i CIA cross section. The 
successive CIA bands corresponding to At; = 0, 1, ...5 are clearly visible in Figures 1 and 2. 
At higher temperatures, weaker transition become more important and fill in the inter-band 
regions. The Av = 1 band centered at ~ 4200 cm -1 (2.4 /jm) is important in T dwarfs as 
it coincides with the K band flux peak, which is caused by a minimum in the combined 
opacities of H 2 and CH 4 . On the other hand, the J and H band flux peaks in the spectra 
of T dwarfs are shaped primarily by H2O and CH4 absorption. Thus, the K/J and K/H 
flux ratios reflect the relative importance of H2 CIA to the H2O and CH4 opacities which is 
affected primarily by the gravity and metallicity of the atmosphere. A relatively weaker K 
band flux indicates a higher gravity (stronger H 2 CIA) or a lower metallicity (weaker H 2 
and CH 4 bands), or both. Our model atmospheres tend to u nderestimate the K band flux, 
resulting in bluer J — K than is observed for late T dwarfs ( ISaumon fc Marleyl 120081 ). We 
have long suspected that this systematic offset comes form an overestimate of the strength 
of the H 2 CIA in the Av = 1 band. 



2.1. New calculation of the H 2 collision-induced absorption 

First-principles calculations of molecular opacities are typically conducted in two steps. 
First, a potential energy surface is calculated for a range of molecular geometries rep- 
resenting variations in bond lengths and bond angles. This potential surface, which is 
mult i- dimensional for polyatomic molecules, is obtained by solving the quantum mechan- 
ical Schrodinger equation for the electrons for each given nuclear geometry using methods 
of quantum chemistry. The potential surface then serves as input for the calculation of 
the eigenstates of the molecule and the induced dipole moment surface from which allowed 
transitions and their oscillator strengths can be calculated to generate a line list. For bands 
originating from excited states, or for relatively high energy transition (in the optical), the 
final state is well above the ground state of the molecule. The calculation of these lines 
requires an accurate knowledge of the potential energy surface far from the minimum that 
approximately corresponds to the ground state energy, which involves many subtle quantum 
mechanical effects and requires extensive calculations for widely separated configurations of 
nuclei. The calculation of higher energy transitions and of "hot bands" is thus increasingly 
difficult. 



The growing evidence from brown dwarf spectra that somet hing was amiss in the At; = 1 
band of H 2 CIA as well as applications to ultracool white dwarfs ( iBergeron. Saumon fc Wesemael 
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19951 ) prompted a new calculation of the H 2 -H 2 and H2-He CIA opacity in the entire regim e 



of temperatures and frequencies of interest (jFrommhold et al.ll2010t lAbel et al.ll201ll |2012| ). 
New, high accuracy potential energy surfaces and induced dipole surfaces were generated for 
these two CIA processes. The resulting opacity was calculated on a fine temperature and fre- 
quency grid, from v = 20 to 10 4 c m" 1 and T = 200- 3000 K for H 2 -H 2 and 20 to 20000 cm" 1 
and T = 200-9900 K for H 2 -He feichard et aDl2012h . The H 2 -H 2 table is being extended to 
the same (T, v) grid as the H 2 -He table. In each case, the entire parameter range is covered 
in a uniform manner, while earlier tabulations used different approximations and potential 
surfaces for different temperature and frequency ranges. As is the case for the older calcu- 
lations, the new H 2 -H 2 op acities are in exce l lent agreement w ith experimental data, which 
are limited to T < 300 K (jFrommholdl Il993l ; lAbel et al.ll201ll ). Comparisons with the CIA 
opacity previously used in modeling brown dwarf atmospheres are shown in Figures Q] and 
[2J As expected, differences are negligible or minimal at low temperatures. The differences 
become substantial at 2000 K in the roto-translational band for H 2 -H 2 (~ 1000 — 2000 cm -1 ) 
and in all roto- vibrational bands (Av > 0). In the fundamental band (Av = 1), which 
affects the K band flux in brown dwarfs, the H 2 -H 2 CIA has decreased by 24% at 1000K 
and by 44% at 2000K, which goes in the direction required to bring the K band flux of the 
models into better agreement with the observations. The blue side of the fundamental band 
of H 2 -H 2 CIA extends into the H photometric band where the absorption has decreased 
somewhat at temperatures above ~ 1500 K. The H 2 -He contribution (Fig. 2) shows similar 
variations, except that the absorption coefficient increases at the peaks of the fundamental 
and first overtone bands. At frequencies above 10 cm" 1 , the corrections can be factors of 3 
or more. This matters less in brown dwarfs as H 2 CIA is not an important source of opacity 
at optical wavelengths and the He abundance is ~ 10% that of H 2 . 

The effect of the change in the CIA opacity is shown in Fig. [3] where spectra computed 
with the old and the new opacity but with the same (T, P) atmospheric structures, are 
compared. All the flux peaks in the near infrared are affected, with rather modest flux 
increases of < 15%. Since the new H 2 CIA opacity is nearly identical to the older calculation 
at T = 300 K and increasingly differs at higher temperatures, the differences in the spectra 
tend to be larger at higher T e g and nearly vanish for T e fj < 500 K except in the Y and J 
bands which are relatively transparent and probe the atmosphere at temperatures well above 
T e ff. Since the CIA opacity increases with the gas density, i.e. with gravity, lower gravity 
spectra are less affected by the new CIA. For example, at logg = 4 (cgs), the flux in the 
peaks of the SED increases by less than 7%. At effective temperatures higher than those 
shown in Fig. [21 a greater change in the spectrum would be expected in cloudless models. 
However, in brown dwarfs atmospheres with T e g > 1400 K the continuum opacity of the 
condensates overwhelms the H 2 CIA opacity and the spectrum is unaffected. 
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3. Ammonia opacity 

In brown dwarf atmospheres, the chemistry of nitrogen is dominated by molecular ni- 
trogen N 2 and ammonia NH 3 , t he latter being favored at higher pressures and lower temper- 
atures ( ILodders fc Fegleyl 120021 ) . The transformation of N 2 into NH 3 occurs at T e g ~ 800 K 
and is, along with the similar conversion of CO to CH 4 , one of the most significant chemical 
changes in brown dwarf atmospheres as T e g decreases. While N2 has no permanent dipole 
moment and is effectively invisible in brown dwarf spe ctral the fundament al band of NH3 
in the 10-11 /xm region is seen in T2 and later dwarfs (ICushing et al.l 120061 ) and the entire 
9-14 jum region is affected by NH 3 absorption ( jSaumon et al.l 120061 ) . Ammonia has a series 
of weaker overtone bands in the near infrared, two of which are centered at 2.0 and 1.5 fim 
and fall in the K and H photometric bands where the SEDs of T dwarfs show emission. 
The apparition of these bands in the spectra of very cool brown dwarfs has been suggested 
as a trigger for a new spectral class (IBurrows et al.ll2003l ; iLeggett et al.ll2007al ; iKirkpatrick 
2008h . 

Saumon et all (boool ) reported the detection of three NH3 features in a R = 3100 spec- 
trum of the T7p Gliese 229B in the 2.02 - 2.045 /im range but other expected NH3 features 
were missing. This tentative detection was followed by a decade where no other near-infrared 
spectrum of T dwarfs was taken at high enough resolution in the K band to reveal NH 3 fea- 
tures. The first unambiguous detection of NH 3 in brown dwarfs was achieved with the 
Infrared Spectrograph of the Spitzer space telescop e that clearly showed the 10-11 Mm band 



in the combined spectrum of e Indi Bab (T1+T6) (iRoellig et al.l 12004 



and then in all dwarfs of type T2 and later observed with the IRS (ICushing et al. 



Mainzer et al. 



7) 



200 



20061 : 



Saumon et al.l 120071 : ICushing et al.l 12008 : IStephens et al.l BOOa ILeggett et al.l l2010bh . Re- 
cently, the search for extremely cool dwarfs has uncovered several objects with T e g < 600 K 
where the blue wing of the H band flux peak appears depressed due the p resence of NH 3 



(iDelorme et al.ll2008t iLucas et al.ll2010l : iLiu et al.ll201ll : ICushing et al.ll201lD . Finally, a new 
R ~ 6000 near-infrared spectrum of the T pff ~ 500 K T9 dwarf UGPS J072227.51-054031.2 
contains NH 3 features ( iBochanski et al.l 120111 ). The detection of NH 3 in the near- infrared 
spectrum of several u ltracool dwarfs is on e of the factors that motivate the creation of the 
new Y spectral class (ICushing et al.l 1201 ll ). 



The analysis of the spectral signature of NH 3 in T dwarfs is complicated by the pres- 
ence of vertical mixing in the atmosphere and by the invisibility of N 2 . Because the time 
scale for the conversion of N 2 to NH 3 is very long at low temperatures, gas eddies from very 



2 Like H2, N2 causes collision- induced absorption ( Borvsow fc Frommholdl 19861 1987 ). but its mole frac- 
tion of ~ 10~ 4 ' 2 makes this contribution to the opacity in brown dwarf atmospheres completely negligible. 
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deep in the atmospheres where N2 and NH3 are in chemical equilibrium at high temper- 
atures can be transported vertically by convection and slow mixing in the radiative zone 
to regions of lower temperatures on a time scale that is much shor ter than the N? — y NH 3 



conversion can take place, resulting in a reduced NH 3 abundance (IFeg 



19961 : lLodders fc Ferievl I2OO2I : ISaumon et aJbood iHubenv fc Burrowsl EoO^ l This process 



ev fc Lodderslll994 



can reduce the observed abundance of NH 3 by factors of 5 — 10, weakening the NH 3 bands 



dwarf observed with the Spitzer IRS with a high enough S/N ratio ( 


Saumon et a . 


2006 


Mainzer et al. 


2007; 


Saumon et al. 


2007; 


Bureasser et al. 


2008b: 


Leggett et al. 


2009, . 


2010b 



and can be fully explained by considering the convective mixing time scale and the time 
scale for N2 — > NH3 conversion. Considering the relatively well understood chemistry of 
nitrogen and the universality of a deep convection zone in the atmospheres of brown dwarfs, 
it is reasonable to assume that the abundance of NH 3 is reduced by vertical transport in all 
T dwarfs. 

In view of the currently very limited opportunities for spectroscopic observations of 
brown dwarfs in the mid-infrared and the relative ease of near-infrared spectroscopy, the 
classification of the growing number of very late T dwarfs and of the new Y dwarfs relies 
primarily on their near infrare d spectral characteristics where NH 3 figures prominently as 



a new spectral index indicator (IDelorme et al. 



Bureasser et al. 


2011) 


Kirkpatrick 


( 


2008). 1 



2008 



Lucas et al 



2010 



Cushing et al. 



Leggett et all ( ]2007ah and 



2011 



indicator of non-equilibrium chemistry, and an important marker for the spectral classifica- 
tion of very late T dwarfs and Y dwarfs demands a high quality line list to compute reliable 
models. 



3.1. New theoretical NH 3 line list 



Our previous models used a NH 3 line list from the HITRAN data base ( IRothman et al. 



2003, 



2005 



2009) supp l emented with room tem perature data in the 6000— 7000 cm 1 range 



( jFreedman et al.l 120081 ; ISharp fc Burrowsl [2007J). In principle, the fundamental band at 
10 — 11 /xm should be well represented in the low temperature (~ 300 K) HITRAN data 
base but the line list becomes increasingly incomplete at shorter wavelengths and higher 
temperatures. The resulting opacity at 400 and 1000 K is shown in Figs. H] and There 
are several gaps in the inter band regions 2.5-2.75 /im, ~ 2.11 /jm and 1.6-1.88 /xm, and 
the list stops at 1.42 /xm. Such a line list is inadequate to analyze the rapidly grow- 
ing spectroscopic data for T e fj < 600 K dwarfs, which has motivated recent experimental 
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2011 



fHa rgreaves. Li fc Bernathll201l!) and theoreti cal work on ammonia (IHuang. Schwenke fc Lee 



Yurchenko. Barber &: Tennyson! 120 111 ) 



A new extensive line list for NH 3 has recently been calculate d from first principles. It 
is based on a new potential energy surface (jYurchenko et al.ll2005l ) that has been refined by 
fitting calculated energy levels to nearl y 400 levels that are accurately known experimentally 
( jYurchenko. Barber fc Tennyson! 1201 ll ). The resulting energy levels agree with the exper- 
i mental values to within 1 cm -1 and generally under 0.1 cm -1 . This 'BYTe' NH 3 line list 
( jYurchenko. Barber fc Tennyson! 120111 ) contains over 1.1 billion lines, compared to ~ 34000 
in HITRAN. It extends from to 12000 cm- 1 (0.83 /an) a nd is applicable up to 1500 K . 
An independent calculation of a NH3 line list is under way ( IHuang. Schwenke fc Leell201ll ). 
Figures H] and show the improvement of the BYTe line list over the HITRAN opacity. At 
400 K, we expect the HITRAN opacity to be quite good since it is based on data acquired 
primarily at room temperature. This is indeed what the BYTe line list shows as both are 
nearly identical near the center of all the bands in common between the two line lists. In 
addition, the BYTe line list includes new bands at 1.65 /zm and shortward of 1.4 fim, and 
contains many weak lines that fill in the inter band regions where there are gaps in the 
HITRAN line list. The differences become more dramatic at higher temperatures (Fig. \5§ 
where the transitions from thermally excited levels come into play, reducing the contrast 
between the peaks and valleys of the opacity and increasing the band opacity by factors 
of ~ 3 — 5 compared to the HITRAN value. The opacity in the 8 — 15 jum also increases 
significantly. 

The effect of the change in NH3 opacity on the synthetic spectra is shown in Fig. [6j 
For this comparison, spectra computed with the HITRAN and the BYTe line lists use the 
same (T, P) atmosphere structure. At T e g = 1200 K (not shown), the spectra are identical 
but at T cS = 800 K there is a small flux decrease of < 20% on the blue side of both the H 
and K peaks. At lower T e g, the chemical equilibrium shifts to strongly favor NH 3 over N 2 as 
the dominant nitrogen bearing species and its opacity becomes more significant compared 
to the competing absorption from H2O and CH4. This leads to dramatic changes in the 
spectrum. All the near infrared flux peaks are reduced by factors of 2 — 4. In particular, 
the new bands at short wavelengths dramatically affect the shape of the Y and J peaks at 
T eff = 400 K. At the extremely low T cff of 200 K, the changes are limited to the spectral 
regions were there are gaps in the HITRAN line list, and lead to very distinctive features 
in the Y and J bands. In particular, the Y band peak is split in two. Figure [6] suggests 
that in addition to the NH 3 — H sp ectral index which m ore or less measures the slope of the 
blue side of the H band flux peak (IDelorme et al.l 120081 ) . another NH 3 -sensitive index could 
be defined to measure the curvature of the J band peak. In the mid-infrared (not shown), 
the most variation is shown by the 800 K spectrum with a < 40% decrease in flux in the 
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9 — 15 fim region. The 400 and 200 K spectra are barely affected in the mid-infrared. All 
these differences decrease somewhat as gravity is increased. The spectra shown in Fig. [6] and 
discussed above are computed in chemical equilibrium. However, we expect that the NH 3 
abundance will be prevented from coming into equilibrium by convection, resulting in surface 
abundances ~ 10 times lower in the upper atmosphere. Consequently, NH3 absorption is 
more muted in non-equilibrium spectra and the changes discussed above are smaller than in 
model that are in chemical equilibrium. The one exception is the T e g = 200 K model which 
is so cold that the entire atmosphere is dominated by NH 3 whose abundance becomes nearly 
insensitive to convective mixing. 



The recently identified Y dwarfs have T eS ~ 300 - 500 K (ICushing et all 120111 ). The 
spectra computed with the new BYTe line list indicate that broad NH 3 features should be 
detectable in the near infrared, most notably as a shoulder on the blue side of the Y band 
peak at T e g ~ 400 K which is split at T c fj- < 300 K by a growing NH3 band centered at 
1.03 (Ltm. The J band peak becomes triangular and the H band peak narrows. The strength 
of these effects is reduced by the non-equilibrium abundance of NH 3 . 



4. Models including both the new H 2 CIA and NH 3 opacities 



We have computed new model atmospheres and synthetic spectra using the new H 2 
CIA and NH 3 opacities. In all o t her aspects, these models are identical to those de- 



scribed inlSaumon &: Marlevl (12 008) ; 



Mar 



ey et al.l ( 120021 ) with additional details provided in 



Cushing et al. ( 20081 ) and" Stephens et al. ( 20091 ) . These new models, which include cloud 



l ess and cloudy atmospheres , wer e used in the r ecent analyzes presented in ICushing et al. 



(120111 ); iLuhman et al.l (120111 ) and lLeggett et al.l (120111 ). These new cloudless model spectra 
are compared to our previous models in Fig. where, in contrast with Figs. [3] and O the 
(T, P) structures are computed with their respective opacities and thus the model structures 
and spectra are self-consistent. In this case, the flux is redistributed along the SED to con- 
serve the total flux (which was not the case in Figs. [3]and[n]). At T c g = 1500 and 1000 K, 
the changes are modest and dominated by the new H2 CIA opacity. The reduced opacity 
in the first overtone band of H2-H2 CIA is responsible for the increased flux in the K band 
peak (~ 15%) and a smaller increase in the H band flux. The T e g = 1500 K spectrum also 
shows a 7% flux increase in the 4 /im peak but the 3-15 fim cloudless spectra are otherwise 
unaffected. Near the T/Y dwarf transition, where T e g ~ 500 K, all four near infrared peaks 
in the spectrum are noticeably affected. Absorption by NH 3 changes the shape of the Y 
and J peaks and shaves off the blue side of the H band peak. Back warming from this NH 3 
absorption pushes the flux out in the K band, which is now ~ 75% brighter, and in the 
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4-8 /j,m region where the flux increases by ~ 30%. 



4.1. Comparison with the spectrum of Gliese 570D 



The T7.5 dwarf Gliese 570D ( iBurgasser et al.ll2000l ) is the most thoroughly studied late 
T dwarf. Optical, near infrared, M band, and Spitzer IRS spectroscopy are available - 
sampling ~ 70% of the total flux emitted - as well as Spitzer IRAC and WISE photome- 
try. As this brown dwarf is a companion to a nearby main sequence star, its distance and 
metallicity are accurately known and its age is fairly well constrained to 2-5 Gyr. Because 
Gl 570D is nearby and relatively bright, the whole data set is of high quality and very 
effectively constrains the models. It is an excellent test object to validate the new mod- 
els and eval uate the improvements brought about by the new opac i ties. Detailed analyzes 
of Gl 570D jGeballe et aDboOll : ISaumon et all bood : beballe et all l2009h have established 
that T eS = 800 — 820 K and log g = 5.09 — 5.23 (cgs), where the ranges of T eS and logg 
are correlated: a higher T eff corresponds to a higher gravity. S olar metallicity was a dopted 
based on the metallicity of the primary ([Fe/H]= 0.09 ± 0.04). ISaumon et al.l (120061 ) estab- 
lished that convective transport drives the NH3 abundance below the value expected from 
chemical equi li brium and obtained an very good fit of the entire SEP of Gl 570D (see also 



Geballe et al.l (120091 )). Figure M shows the spectrum computed by ISaumon et al.l (120061 ) for 
those parameters within the above range that best reproduce the data and a new model 
computed with the same atmospheric parameters but with the new H2 CIA and NH3 opac- 
ities. There is significant improvement in the K band peak which is now fitted very well. 
There is also a modest improvement in the 5.5-7.5 fim region and beyond 9 fim where NH 3 
absorption features overlap with H 2 absorption. The new H 2 CIA calculation, which is pri- 
marily responsible for the increased flux in the K band represents a significant improvement 
over the previously available tabulation. Outside of the 9-14 /iin region, NH3 absorption is 
rather weak at this T e g, and the new NH3 opacity leads to very modest changes. We believe 
that the remaining discrepancies in our fot of Gl 570D, such as the width of the J band 
peak and the depth of the 1.6 /zm band of CH 4 , are primarily due to the very incomplete 
methane line list, which is prese ntly in a similar state as the ammonia line list was prior to 
the first-principles calculation of lYurchenko. Barber fc Tennyson! (1201 11 ). 



4.2. Comparison with near-infrared photometry 



While a comparison with the SED of Gl 570D provides a detailed test of the new 
models, the overall effect of the new opacities can be better appreciated in color-magnitude 
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diagrams where synthetic colors can be compare d to a large set of photom etric data. Prior 
to 2003, our models used the solar abundances of lLodders fc Fegleyl (119981 ) and the cloudless 
sequence was in gener ally good agreement with the near infrared colors of late T dwarfs 
( Burgasser et al. 2002^. A significant downw ard revision of the CNO abundances in the 
solar photosph ere (lAllende Prieto et al.ll2002l ) motivated a change of the solar abundances 
in our models ( Lodders 20031 ) which moved the cloudless sequence to the blue in J — K. The 
resulting sequence of [M/H]=0 cloudless models are too blue compared to late field T dwarfs, 
particularly in J — K, when compa r ed to the near infrared colors of late-T field dwarfs (Fig. M 
and Fig. lOe of lSaumon fc Marlevi (bo08h )R The new H2-H2 CIA has a weaker first overtone 



band which increases the K band flux and moves the cloudless models to redder J—K, largely 
compensating for the effect of the reduced CNO abundances (Fig. [9]). The new sequence of 
models is still slightly bluer by < 0.15 than the late T dwarf sequence which may be due 
to missing CH 4 opacity on the J band peak. The spread of the observed sequence of late T 
dwarf colors is significant an d can be accounted f or by the observed distribution in metallicity 
of field dwarfs (Fig. 15f of ISaumon fc Marleyl (120081 ) ). The very latest T dwarfs recently 
discovered, with Mk ^ 17 reveal a saturation in J — K and are redder than the calculated 
cloudless sequences. This may be caused by several limitations in our models such as the 
near infrared opacity of CH 4 , the presence of clouds of compounds with low condensation 
temperature such as Na2S, or the appearance of water clouds. Note that water condensation 
is included in the chemistry of the atmosphere and the gas phase H2O abundance is reduced 
in the low-T e g models where water condenses. We are investigating the role of clouds in 
low-T e fj atmospheres and the results will be reported in a separate publication. On the other 
hand, our cloudy sequences account for iron a nd silicate clouds and reproduce the SEDs and 



infrared colors of fie ld L dwarfs fairly well ( ISaumon fc Marleyl 120081 ; iCushing et al.l 12008 



3 p: 
8; 



Stephens et al.l 120091 ) . Cl oudy models with T e ff > 900 K are representative of field L dwarfs 
and the HR 8799 planets (IMarley et al.ll2012l ) and are shown in Fig. |9j These hotter, cloudy 
models are barely affected by the new opacities and their near infrared colors are very nearly 
identical to the older ones down to T c g = 900 K. This is a consequence of the opacity of 
condensates which reduces the importance of the H2 CIA and of the higher temperatures of 
the cloudy models which reduce the abundance of NH3. 

The same general picture emerges from the Mj versus J — H color-magnitude diagram 
(Fig. [10]) although the new cloudless model now represents the J — H sequence of the late 
field T dwarfs quite welfl The data suggest that the J — H color saturates around —0.6 (see 



3 Note that the clou dless COND models of lAllard et al.l (|2001l ) use the older solar abundances of 



Grevesse fe; Noelsl (|1993l ) and agree very well with the observed near-infrared colors of late T dwarfs. 



4 This may be somewhat fortuitous in view of the model limitations mentioned above, particularly the 
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also iKirkpatrick et al.l ( 1201 ll ) at Mj ~ 17 which is further supported by the position of the 
Y0 dwarf W ISEP J1541 - 2250 in t his color-magnitude diag ram (Mj = 23.9 ± 0.8, J-H = 
0.17 ± 0.63, ICushing et al.l (1201 lh : IKirkpatrick et al.l ( 1201 ll )). Our cloudless sequences do 
show such a minimum in J — H at T e fj ~ 400 K and a rapid turnover to the red at lower 
T e ff. However, the models do not turn over to the red soon enough and reach values of 
J — H ~ —1, which is bluer than is observed. It will be interesting to see how the near 
infrared color-magnitude diagrams will be affected by a more complete CH 4 line list when it 
becomes available, or by the inclusion of a water cloud in the models. 



5. Ammonia and Y dwarfs 



Increased absorption in the near infrared has been observed on the blue side of the H 
band peak in the very late T dwarfs and in the new Y dwarfs and it is thought to be due 
to NH 3 . The str ength of the 1.5 /xm b and of ammonia can be quantified with the NH 3 — H 
index defined by iDelorme et al.l ( 120081 ) 



NFT — H 



(2) 



Smaller values of NH3 — H indicate a relatively low flux on the blue side of the H band 
peak that is expected from stronger NH3 absorption. A compilation of the data for the 
latest dwarfs known is shown in Fig. [TTJ The NH 3 — H index strong l y cor relates with 
the width of t he J band peak (measured by the Wj index, I Warren et al.l (120071 )). and with 
spectral type (IBurningham et al.ll2008l ; ICushing et al.ll201ll ). with later types having smaller 
NH3 — H. Figure 1TT1 sho ws a break in the slope at Wj ~ 0.3 corresponding to spectral type 
T8 ( Cushing et al.l 2011 ). Based on the appearance of the blue side of the H band peak of 
the Y0 dwarf WISEP J1738+ 2732 - one of the m ost extreme objects in Fig. HU- compared 



to earlier T8 and T9 dwarfs. ICushing et al.l (120111 ) tentatively ascribe the absorption to NH3. 



The models shown in Fig. [TT] reproduce the observed trend very well and indicate that 
the reduced abundance of NH3 caused by non-equilibrium chemistry is favored. The dashed 
lines show the indices computed from the same models after the opacity of NH3 has been 
removed with a trend that is strongly at odds with the data. The decrease in the modeled 
NH 3 — H at T cS < 600 K (or spectral type later than T8) is caused by the 1.5 /xm band of 
NH 3 . The expectation that NH 3 bands would appear in the near infrared at low enough T e g 
and its detection as a depression in the blue side of the H band peak is thus well supported 



poorly modeled 1.6 /im band of CH4 which affects the H band (Fig. [5]) 
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on a quantitative basis by the models. However, the non-equilibrium models still show a 
systematic offset from the data, even for the low gravity of log g = 4 (cgs) which, for these 
low T cff , would correspond to a rather low mass of ~ 10 Mj upiter . It is very likely that the 
incomplete CH 4 line list is responsible for this systematic shift since methane is a strong 
absorber in the H band and a moderate absorber in the J band (Fig. [8]). 

Without the benefit of models, or with models of limited fidelity, the interpretation 
of the observed trends in spectral indices can be ambiguous, particularly when considering 
their variations over a wide range of T e g. On the other hand, the direct detection of spectral 
features provides much more secure evidence of the presence of specific molecular species. 
While the detection of NH 3 in the 2.02-2 .045 um region o f the s pectrum of the T7p Gl 
229B (jSaumon et all 120001 ) was tentative, iBochanski et all (1201 ll ) identified 11 ammonia 
features in their medium resolution (R ~ 6000) near infrared spectrum of the T9 dwarf 
UGPS 072227.51-054031.2 (here after UGPS 0722-05). Figure QJ shows the portion of 
their spectrum with the most NH3 features. For comparison, the top curve shows the log of 
the NH3 opacity computed from the BYTe line list at the conditions corresponding to the 
atmospheric level where this part of the spectrum is emitted. The lowest curve is an unsealed 
model spectrum wit h T eff = 500 K and log g = 4.25, solar metallicit y, and K zz = 10 4 cm 2 /s 
( ILeggett et al.ll201ll ). The six NH3 absorption features identified by IBochanski et all (1201 ll ) 
in this part of the spectrum all have corresponding troughs in the model spectrum. Four of 
these features also match peaks in the NH 3 opacity but the other two (1.5264 and 1.5411 /zm) 
do not and their attribution to NH 3 is not secure. 

To better identify NH3 absorption features in the spectrum of UGPS 0722—05 we show in 
Fig. [12] a spectrum computed from the same atmosphere model but without any NH 3 opacity 
(blue curve). There is obviously much NH 3 absorption over the entire range shown, even 
considering the fact that the model spectrum has a reduced (non-equilibrium) abundance 
of NH 3 . This is fully consistent with the smaller value of the NH 3 — H index of UGPS 
0722—05 compared to earlier T dwarfs (Fig. [TT]) . However, a careful comparison of the 
spectra with and without NH 3 does not reveal any distinct feature at this spectral resolution. 
Every feature in one spectrum appears in the other, albeit with a different amplitude. For 
example, the 1.5142 /im absorption feature matches a peak in the NH 3 opacity and an edge 
in the model spectrum (Fig. 112"]) . but that edge also appears in the spectrum computed 
without NH 3 , making the attribution doubtful. The model spectrum reproduces individual 
features of the observed spectrum fairly well (after scaling), but not exactly, which makes 
the attributio n of specific feat u res p roblematic. We find this to be also the case in other 
regions where IBochanski et all ( 1201 ll ) found NH 3 features (1.21— 1.27 jt/m, 1.56-1.64 /zm and 
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at 1.990 /mi) @ 

The situation is different in the 1.96-2.05 /im spectral region where ISaumon et al.l (120001 ) 
found 3 ammonia features in Gl 229B. Figure [TBI shows this spectral region with the same 
models and opacity shown in Fig. [12] (with different scalings for clarity). The vertical solid 
and dashed lines correspond to 19 NH3 absorption features that are clearly identifiable from 
the difference between the spectrum computed with the non-equilibrium NH3 abundance 
(red) and the spectrum computed without NH 3 opacity (blue). In this spectral range, almost 
every peak in the NH 3 opacity (green) corresponds to a distinct absorption feature in the 
synthetic spectrum. The solid lines identify 8 absorption features that are present in the 
spectrum of UGPS 0722-05§, and the dashed lines indicate 6 NH 3 features that are not 
clearly detected in the spectrum and 5 that appear to be missi ng. Two of the detected 



features (2.0374 and 2.0455 /mi) match those found in Gl 229B by ISaumon et al.l ( 120001 ). An 
absorption feature clearly identified in iBochanski et al.l ( 1201 it ) at 1.9905 /mi (dotted line) 
corresponds to a well-defined absorption feature in the model spectrum. This feature is 
barely affected by the removal of all NH3 opacity, however, and does not match very well 
with the NH 3 opacity peak centered at 1.9896 /im. In our model spectrum, this absorption 
is caused by H 2 0. 

We conclude that the detection of ammonia absorption features in the near infrared 
spectra of very late T dwarfs and Y dwarfs is easier in the K band near 2 /mi than in the 
H band near 1.53 /im. The main reason for this is that the NH 3 opacity in the K band 
is composed of a series of more or less equally spaced and well separated peaks (Fig. [T3"|) 
while in the H band it has a denser irregular structure with a smaller dynamic range as 
well as a peak opacity that is about one order of magnitude smaller (Figs. [12] and H]). With 
higher resolving power it becomes possible to isolate NH 3 features in the H band. For 
example, with R ~ 20000, about a dozen features stand out in the 1.51— 1.57 /im region of 
the model shown in Fig. IT2"1 Fo r dwarfs cooler than ~ 400 K, the condensation of water 



becomes significant ( IBurrows et al.ll2003l ). removing H2O from the gas and unmasking the 
NH 3 bands. For instance, the blue side of the H band peak becomes completely dominated 
by NH 3 absorption and NH 3 features should be readily detectable in the spectrum. This is 
also true in the 2 - 2.05 /im region, although the flux is much lower than at 1.55 /im. The 



5 For completeness, we note that the t wo features in the low-resolution (R ~ 500) spectrum of UGPS 
0722—05 pointed out bv lLucas et al.l (|2010r ) at ~ 1.275 /^m and at ~ 1.282 /im can be identified with synthetic 
spectra at higher resolution as blends of H2O features. A weaker blend of H2S features also contributes to 
the latter. 



6 If we ignore the single high- valued pixel making up the spike in the data at 2.0175 /im, the NH3 feature 
stands out. 
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NH3 bands will thus become stronge r with lower T c g until ammonia itself begins to condense 



at T eS ~ 160 K f lBurrows et all 12003 ) . 



6. Conclusion 

We have computed new cloudy and cloudless models and synthetic spectra0 of brown 
dwarf atmospheres using recent first-principles calculations of the opacity of NH3 and of the 
collision-induced absorption of H 2 . Both constitute very significant theoretical improvements 
in the modeling of the opacity and, in the case of NH 3 , the BYTe line list greatly expands 
over the HITRAN line list that is based primarily on experimental results. The cloudy 
models, relevant to L dwarfs, are barely affected by the new opacity because of their higher 
T e Q where NH 3 has a very low abundance and of the condensate opacity that largely conceals 
the H 2 CIA @. We have thus focused on the role of these two opacities in cloudless models, 
which are relevant to spectral types ~ T4 and later. 

The most notable change in the collision-induced absorption is a reduction of the opacity 
in the first roto- vibrational band centered at 4200 cm -1 (2.4 fim), resulting in an increased K 



band flux. This corrects a deficiency in our earlier fit of the SED of Gl 570D (jSaumon et al. 



20061 ) and we now obtain a remarkable fit of nearly the entire SED of this late T dwarf (Fig. 
[E]). Another consequence is that the cloudless sequence shifts toward bluer colors in near 
infrared color-magnitude diagrams and is in much better agreement with the colors of the 
late T dwarfs in the field. This new accurate H2 CIA opacity will strengthen the analysis of 
gravity and metallicity indicators in the spectra of late T dwarfs and L subdwarfs and help 
reduce the uncomfortably large uncertainty in spectroscopic determinations of the gravity of 
brown dwarfs. 

The new BYTe line list for NH3 greatly improves the near infrared opacity of ammonia 
by adding several new bands at short wavelengths as well as a very large number of "hot 
transitions" that fill in gaps in the HITRAN line list. We can now reliably calculate the 
opacity of NH 3 at temperatures of several hundreds of degrees, which is crucial even for the 
coolest known dwarfs as the Y, J and H bands typically probe temperatures well above 
T e ff. The importance of this line list for the study of the T/Y transition and the new Y 
dwarfs cannot be overstated. Our models show that NH 3 develops very distinct signatures 
in the Y and J band flux peaks, and confirms the tentative identification of NH 3 absorption 



7 The synthetic spectra are available upon request to D. Saumon (dsaumon@lanl.gov). 



"This is no longer the case in low-metallicity L subdwarfs where H 2 CIA is noticeable ([Burgasser et al 
2003al l2009t ICushing et al~ll2009h . 
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on the blue side of the H band peak which had been anticipated by iBurrows et all ( 120031 ). 
More quantitatively, the new models reproduce the trend in the spectral indices NH 3 — H 
and Wj, confirming that the rapid decrease in NH 3 — H at spectral types later than T8 is 
due to NH 3 absorpti on. Thus, all the new WISE objects assigned to the Y dwarf class by 
Cushing et all (120111 ) show NH3 absorption in the H band. Furthermore, the data favors 
models where vertical transport reduces the NH3 abundance in the observable part of the 
atmosphere. At a more detailed level, the medium resolution near infrared spectrum of the 
T9 UGPS 0722-05 (IBochanski et al.ll201ll ) provides a rare opportunity to compare data with 
the detailed predictions of spectral features by the models and of t he underlying mo l ecular 
opacities. While we found that none of the features pointed out by IBochanski et all (1201 if ) 
can be unambiguously attributed to NH3, our new models allow the identification of 9 other 
NH 3 absorption features around 2 /im. The models clearly show that while there is strong 
NH 3 absorption in the H band in late T dwarfs and in the Y0 dwarfs, the general behavior 
of the NH 3 absorption in the 1.55 /zm region makes it more difficult to isolate individual 
features than in the 2.0 fim region, requiring a resolving power of the order of R ~ 20000 
compared to ~ 6000 in the K band. 

Despite the progress that we have demonstrated here, obvious problems remain with 
the models. For instance, the fit of Gl 570D still has discrepancies that stand out. The 
near infrared color-magnitude diagrams for field dwarfs fail to reproduce the saturation (and 
possible turnover) revealed by the data at very low T c g. The spectral indices NH 3 — H and 
Wj do not yet provide a quantitative diagnostic of the atmospheric structure and chemistry. 
With these two new calculations of molecular opacities that affect primarily late T dwarfs 
and the new class of Y dwarfs, reliable opacities are available for all important gas phase 
absorbers except for CH4. The inadequacy of the the current line list o f methane has been 



much maligned ev er since the very first spectral modeling of a T dwarf ( iMarley et al.lll996 



Allard et al.lll996l ). We believe that for cloudless atmosphere models, the CH 4 opacity is the 
only remaining obstacle of significance to modeling the spectroscopic data accurately. The 
availability of a CH4 line list comparable in quality and scope to the BYTe NH 3 line list 
would open the door to reliable determinations of T dwarf gravities and metallicities, two 
fundamental parameters that have remained elusive for the vast majority of brown dwarfs. 
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Fig. 1. — Collision-induced absorption coefficient for H 2 -H 2 collisions at T = 500, 1000 and 
2000 K, from bottom to top, respectively. Each broad peak corresponds to a change in the 
vibrational quantum number of Av = 0, 1, 2... from left t o right, res p ective ly. The calcula- 
tions of Borysow and collaborators (jBorysow et al.l (120011): iBorvsow d2002l), and references 



there in) and the recent work of Abel & Frommhold ( IFrommhold et al 



2010: Abel et al. 



201 ll ) are shown in red and blue, respectively. 1 amagat= 2.6867774 x 10 19 cm 3 . [See the 



electronic edition of the Journal for a color version of this figure.} 
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Fig. 2. — Same as Figure [ D for H?-He CIA . The older calculation labeled "Borysow" is 
from Borysow et al. ( 1989 ). Ijorgensen et al. J 2000 ). and references therein. The "Abel & 
Frommhold" calculation is from Abel et al. ( 20121 ). [See the electronic edition of the Journal 
for a color version of this figure] 
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Fig. 3. — Effect of the new H2-H2 and H\2-He CIA opacity on synthetic spectra of brown 
dwarfs. The spectra shown are cloudless models with T e g = 1500, 1000 and 500 K, with 
log g = 5 (cgs) and solar metallicity. The spectra computed with the new CIA opacities 
are shown in blue. The red lines show spectra computed with the older CIA opacity and 
the same (T, P) structures. The fluxes are calculated for d — 10 pc and are displayed at a 
resolving power R = 500. 
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Fig. 4. — Absorption coefficient of NH3 at T = 400 K and P = 1 bar. The red line shows the 
opacity calculated from the line list that we have used in our previous atmosphere models. 
It is primarily based on the HITRAN data base supplemented with more recent laboratory 
measurements (IFreedman et al.l 120081 ). The blue line labele d "BYTe" shows the absorption 
coeffic ient computed from the new first-principles line list of lYurchenko. Barber fc Tennyson 
(1201 ll ) . For clarity, the absorption coefficient is shown at a resolving power of R = 200. 
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Fig. 5. — Same as Figure H] for T = 1000 K. 
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Fig. 6. — Effect of the new NH3 opacity on near infrared spectra of cool brown dwarfs. Three 
models are shown with T e g = 800, 400 and 200 K. All models shown have solar composition, 
log g = 4, are cloudless and in chemical equilibrium. Models computed with the HITRAN 
line list for NH 3 are shown in red while those computed from the same (T, P) structure but 
with the BYTe NH3 line list are in blue. The spectra are shown as absolute fluxes for a 
distance of 10 pc and at R = 500 for clarity. 
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Fig. 7. — Comparison of spectra including the new H 2 CIA and the BYTe NH 3 opacities 
(blue) with spectra from a previous generation that only differ in these two sources of opacity 
(red). These cloudless spectra have logg = 5 (cgs), solar metallicity and are computed in 
chemical equilibrium. The effective temperatures are 1500, 1000 and 500 K. In contrast with 
Figs. |3] and El the atmospheric (T, P) structures have been computed with the appropriate 
(old or new) opacities. The spectra are plotted at R = 500 for clarity. 
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Fig. 8. — Spectrum of Gliese 570 D. The data are shown in black (IBurgasser et al.l 12003b 



Geballe et al.ll200ll ; ICushing et al.l 12006 



Geballe et al.l 120091) . The spectrum obtained with 



our previous models is shown in red (jSaumon et all2006l ). The blue curve shows the spectrum 
computed with the new H2 CIA and the BYTe NH3 line list with the same model parameters: 
T eff = 820 K, logg = 5.23 (cgs), [M/H]=0, and eddy diffusion coefficient \ogK zz (cm 2 /s)=4.5. 
The models are plotted at a resolving power of R = 500 (top panel) and R = 100 (bottom 
panel), approximating that of the data. The models are not normalized to the data but are 
fluxes at Earth computed with the known distance of Gl 570D and the radius obtained from 
the evolution with the above parameters. 
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Fig. 9. — Mx versus J — K color-color diagram (MKO system) comp aring models with th e 
near infrared photo metry of field MLT dwarfs. The photometry is from lLeggett et al.l ( 120021 ). 
Knapp et al.l (120041 ). and iMarocco et al.l ( 120101 ) with M dwarfs in black, L dwarfs in red and 
T dwarfs in blue. All known binaries have been removed from the sample except those with 
resolved MKO photometry: e Indi B JMcCaughrean et~al~1l2004h . SDSS J102109.69-030420 1 



and SDSS J042348.57-0 41403.5 (Bureasser et alil2006l). and Kelu-1 dLiu fc Leeeett 



The parallaxes ar e from IPerrvman et al. ( 1997 ). Dahn et al. ( 2002 ). Tinnev et al. 



Vrba et al.l ( 120041 ). IMarocco et al.l ( 120101 ) . and various references in iLeggett et al 



2005). 



(2003), 



(12002( 1 . 



Solid blue curve: Cloudless models calculated with the new H2 CIA and the BYTe NH3 line 
list and an eddy diffusion coefficient K zz = 10 4 cm 2 /s (e.g. these model are not in chemical 
equilibrium). The curve extends from T e g = 1500 K to 500 K, with open circles indicating 
T e ff = 1500 K and 1000 K. The dashed blue curve shows m odels computed with old er H 2 
CIA and NH3 opacities and K zz = (chemical equilibrium) dSaumon fc Marlevl 120081 ) . The 



dashed red curves are cloudy models in chemical equilibrium from ISaumon fc Marleyl (120081 ) 
extending from 2400 K to 900 K with open circles indicating (from the top) T e g = 2000, 
1500, and 1000 K. Each curve corresponds to a different value of the cloud sedimentation 
parameter (see fc MarU » Lemon fc Marleyl «) of / sed = 1, 2 and 3, 

from right to left, respectively. Cloudy models computed with the new opacities overlap the 
cloudy sequences shown. All models shown have solar metallicity and log g = 5. 



-34- 




Fig. 10.— Same as Figure M for Mj versu s J - H. The Y0 dwarf WISEP J1541-2250 with 
Mj = 23.9±0.8 and J-H = 0.17±0.63 jCushing et aJboiltEcirkpatrick et alJboilh . falls 
well outside this figure. 
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Fig. 11. — Spectral indices showing the growth of absorption in the blue wing of the H 
band of lat e T dwarfs. T h e spe ctral indices NH 3 — H and Wj are defined in iDelorme et al. 



( 120081 ) and IWarren et al.l (120071 ). respectively. The curves show the indices computed from 
cloudless synthetic spectra starting at T e s = 200 (lower left corner) to 1000 K (right) with 
solid dots marking steps of 200 K. Two gravities, log g = 4 (cgs, red curves) and 5 (blue), are 
shown. Dotted lines show models in chemical equilibrium and solid curves are models driven 
out of equilibrium by convective transport and mixing in the radiative zone parametrized 
by the eddy diffusion coefficient K zz = 10 4 cm 2 /s. The dashed lines show spectra computed 
with the same ( T, P) structures after removin g the NH3 opac i ty. T he dwarfs Gl 570D (T8, 



Delorme et al.l (12008^. U GPS 0722-05 (T9, ICushing et al.l tOllh). W ISEP J1738+2732 



(Y0, ICushing et al.l (1201 lh. and W 
Data compi led from 



2009 



Warren et al. 



SEP J 1 541-2250 (Y0, 
fioj; 



Delorme et al. 



Cushing et al.l (|2011l0 are labeled . 



2011al Jbh: ICushing et al.l (boilh . and Liu et al.l J201lL 



(120081 ): iBurningham et al.l (120081 . 
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Fig. 12.— Medium resolution spectrum (R ~ 6000) of UGPS0722-05 (iBochanski et al. 



20 111 ) (black). Th e observed spectrum ha s been rescaled, shifted vertically, and the radial 
velocity (47km/s, iBochanski et al.l (120111 )) has been removed. A cl oudless synthetic sp ec- 
trum, with T eS = 500 K, logg = 4.25, [M/H]=0, K zz = 10 4 cm 2 /s ( TLeggett et aDboilh is 
shown in red. The blue curve represents the same spectrum calculated without any NH 3 
opacity. The log 10 of the NH 3 opacity at 650 K and 1 bar, corresponding to the emission 
level of this part of the spectrum in the model, is shown in green at an arbitrary scale. The 
model spectra and the opacity are shown at the spectral resol ution of the data . The vertical 
dotted lines show the NH 3 absorption features identified by IBochanski et al.l (120111 ) in the 
spectrum of UGPS 0722—05 which should correspond to peaks in the NH 3 opacity. 
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Fig. 13. — Same as Fig. [12] but in the K band. The models have been rescal e d for clarity. 
The vertical dotted line shows the NH3 feature identified by iBochanski et al.l ( 120111 ) in the 
spectrum of UGPS 0722—05. The other vertical lines are NH 3 absorption features predicted 
by the model spectrum, which all match peaks in the NH 3 opacity. Those that are also 
apparent in the data are shown with solid lines while those that are ambiguous or missing 
are shown with dashed lines. 



